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I.   INTRODUCTION 

As  weapons  systems  grow  in  complexity,  the  cost  in  time  and 
money  of  developing  and  testing  advanced  weapons  is  so  great  that  mili<- 
tary  planners  must  insure  that  testing  procedures  for  performance  and 
tactics  are  formulated  so  as  to  obtain  reliable  data  for  each  test. 
Thus  the  development  of  valid  and  reliable  models  for  test  design  is  a 
matter  of  great  importance.   Model  manipulation  is  a  far  cheaper  method 
for  developing  good  system  design  than  actual  trial-and-error  tests. 
This  paper  presents  a  stochastic  approach  to  a  combat  model  of  tanks  vs. 
tanks  with  emphasize  on  the  following: 

1.  A  scheme  for  development  of  trial  validation  criteria. 

2.  A  useful  procedure  that  can  be  expanded  to  a  Semi- 
Markov  model  for  combat  (land,  sea  or  air  combat). 

3.  Approachs  to  data  presentation. 

4.  Insight  into  data  analysis  and  data  reduction. 

We  let  a  point  P,  with  coordinates  (k,£)  represent  the  number 
of  casulties  of  two  opposing  forces.   Thus  the  coordinates  are  integers 
0  <  k  <  N   and  0  <  £  <  M,  where  N  is  the  size  (number  of  units)  of 
the  friendly  force  and  M  is  the  size  of  the  threat  force.   It  is  con- 
venient to  normalize  this  rectangle  to  the  unit  square,  letting  x  =  k/n, 
y  =  £/m  (see  figure  1).   The  trial  events  are  thus  points  on  the  unit 
square.   Given  a  point  (x,y)  for  any  given  weapon  system  there  is  some 
probability  of  transition  to  some  other  point  (x  ,y  ),  as  the  result  of 

combat.   The  trial  (battle)  is  over  when  either  force  becomes  too  small, 
that  is,  our  point  reaches  some  predetermined  barrier  on  the  top  or 
on  the  right  side  of  the  square.   It  is  our  belief  that  if,  in  any  real 
experiment  (mock  battle),  one  force  quickly  overwhelms  the  other,  the 
experiment  is  a  failure,  in  that  it  did  not  really  test  the  system. 
Graphically  the  point  (x,y)  would  move  nearly  vertically  or  nearly 
horizontally  across  the  unit  square.   Thus  we  wish  to  study  how  to 
design  the  experiment  so  that  the  point  (x,y)  will,  with  some  reason- 
able certainty,  remain  in  a  certain  size  cone  about  the  45°  diagonal. 

In  describing  our  procedure,  semi-Markov  process  and  random 
walk  terminology  are  used  throughout.   We  are  concerned  with  the  develop- 
ment of  the  model  and  its  ability  to  present  data.   In  particular  we 
will  make  certain  simplifying  assumptions  concerning  the  transition 
probability  matrix  for  illustrative  purposes. 

The  technical  portion  of  this  paper  seems  new  and  provides  a 
-discrete  state  two-dimensional  stochastic  model- of  land,  sea  or  air 
combat. 
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II  THE  MODEL 

We  assume  that  the  duration  of  a  trial  is  a  random  variable 
T  and  AT  denotes  the  time  interval  used  to  update  changes  in  systems' 
state  (friendly  and/or  enemy  kills).   The  realization  of  times  where 
changes  in  system  states  are  being  updated  are  expressed  in  terms  of: 

0,  AT,  2(AT),  3(AT),  4 (AT),  ... 

AT  will  be  selected  based  on  the  following  factors: 

1.  AT  cannot  be  selected  too  "small"  due  to  limita- 
tions of  the  instrument  used. 

2.  AT  cannot  be  selected  so  "large"  as  to  permit  more 
than  one  system  to  be  attrited  from  each  force. 

After  AT  is  selected,  the  particle  (representing  the  systems' 
state)  may  change  its  position  or  stay  put  as  time  progresses  depending 
upon  the  results  of  some  corresponding  exchange  of  fire,  but  once  it 
reaches  a  specified  point  on  the  lines  x  =  b  and/or  y  =  b   (where 
50%  <  b  <  70%),  the  system  or  systems  are  considered  to  be  saturated 
and  the  problem  (engagement-trial)  .therefore  terminated.   This  termina- 
tion process  compares  with  actual  battle  conditions  where  a  specified 
force  is  considered  noneffective  when  a  percentage  (50-70%)  of  casualties 
are  sustained.   (See  Figure  1) 
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FIGURE  1, 


x-axis 


x4- 


4 


Casualties  (%)  of  Friendly  Force 

At  the  commencement  of  the  engagement  (trial),  time   t,   the 
system  (friendly  and  aggressor)  is  in  state  (x,y).   At  the  conclusion 
of  a  AT  interval  the  position  of  the  particle,  say  (x,y)  moves  to  a 
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position  denoted  by  S  ,  S  ,  and  S    as  outcomes  of  the  AT  time 
y  y   x      xy 

interval  when  the  particle  moves  from  (x,y)  to         l\ , 

X'  y   M 


x  +  h  y  +  m 


*  +  j.  y  ' 


respectively  (see  Figure  2).   S   can  be  used  to  further 


denote  the  status  of  the  particle;<given  no  change  in  the  state  or  incre- 
ment of  the  particle.   Based  on  the  hypotheses  used  in  this  model,  all  the 
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FIGURE  2. 


possible  positions  of  the  particle  are  in  the  unit  square  and  no 

outcomes  are  possible  other  than  the  above  four,  i.e.,  S  ,  S  ,  S  ,  S 

v  o   x   y   xy 

However,  instead  of  using  the  unit  square  as  the  sample  space,  it  is 

easier  to  return  to  the  rectangle  fi: 

fl  -  {(x,  y);  x  =  0,1,2,. ..,N;   y  =  0,1,2, ... ,M} 

Further,  we  define  the  sets: 

A  -  {(x,y) :   0  <  x  <  N,  y  =  M} 

B  -  {(x,y):   x  =  N,  0  <  y  <  M} 

C  =  {(N,M)} 

1=  {(x,y):  0<x<N,  0<y<M} 
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We  shall  call  the  totality  of  points  in  A  or  in  B  or  in  C  the 
absorbing  barriers,  while  all  the  points  in  I  are  called  interior 
points.  ( 

For  each  point  (x,y)  in  I,  we  assume  we  are  given  the 
probabilities: 

Pr{Sy}  =  p,  p  >  0 

W  =  q'  q  >  ° 

P  {S   }=  r  and  P{S}=s,  r,  s>0 
r  xy         r  o 

such  that  p+q+r+s=l  and  they  are  independent  of  T.   For  sim- 
plicity assume  also  they  are  independent  of  (x,y).   In  the  actual  experi- 
ment, these  probabilities  are  functions  of  location  (x,y).   That  is  to 
say,  the  transition  probabilities  depend  upon  the  remaining  forces  left 
in  the  battle  zone.   If  one  is  interested  in  developing  a  stochastic 
model  for  a  combat  experiment,  the  probabilities  p,  q,  r,   and  s  will 
generally  be  functions  of  time,  locations  of  the  particle,  and  remaining 
units  left  in  the  battle  zone. 

For  each  point  in  the  absorbing  barriers,  we  assume  that: 

P  (S  }  =  1, 
r  o 

P  {s  }  =  P  {S  }  =  P  {S  }  =0. 
r  x      r  y      r  xy 

Denote  P(x,y|t)  the  probability  that  the  particle  for  the  first  time 
reaches  the  point  (x,y)  at  time  t  and  denote  P(x,y|t)   the  probability 
that  the  particle  is  in  state  (x,y)  at  time  t,   and  let  i(x,y)   be 
the  first  passage- time  through  the  point  (x,y)  namely,  the  number  of 
(AT)'s  to  be  performed  until  the  particle  for  the  first  time  reaches  the 
point  (x,y).   Then,  for  any  assigned  values  of  N  and  M,   i(x,y)  may 
be  considered  as  a  random  variable,  the  realization  of  which  equals 
one  of  the  following  values,  i.e., 


where: 


Denote: 


•t,t  +1,  t  +2,  ...  ,  t  +  v,...  (0  %   v  <  ») 

o'  o      o         '  o 


t  =  max  (x,y)  . 


K  (x,y,t)  =  {k;  max(x  +  y  -  t,  0)  <,   k  <  min(x,y)} 
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Denoting  by  P  (i,j,k,£)   the  probability  function  of  the  multinomial 

t 
distribution  (p  +  q  +  r  +  s)  ,  the  following  well-known  identity  holds: 

Pt(i,j,k,JO  =  1;j;fciAt  P1  ^  ^  sl'>    (i+J+k+^  =  t) 
and  to  any  interior  point  (x,  y) ,  i.e.,  (x,y)cl  , 

P(x,y|t)  =  I       t\ x-k  y-k  k   t+k-x-y 

keK  (s-k)! (y-k) Ik! (t+k-x-y) I   q    p    r  S 

for   t  =  t  ,  t  +  1,  t  +2,  ...,  where  t  =  max(x,y). 
o   o       o  o 

For  (x,y)  and  t  such  that   t  <  max(x,y),  it  is  clear  that: 

P(x,y|t)  -  0 

For  the  first-passage  probability  p(x»y|t)  we  have  the  follow- 
ing equations: 

P(0,0|0)  -  1  , 

P(x,y|t)  =  P(x,y|t)  -  sP(x,y|t-l)  . 

If  (x,y)el,  we  should  have: 

■o/        1^    V  (x+y-k) !      x-k  y-k  k 

p(x,y|t)  =  J  (x-k): (y-k):  k.'  q      p      r 

keK(x,y,t)       w 


f   t-1  I 
x+y-k-1 


t+k-x-y 

s      ' 


provided  t  >  t  . 


Proofs  of  all  formulas  and  optimal  mixes  of  weapon  systems  used 
for  both  the  friendly  force  and  the  threat  force  will  be  treated  in  a 
separate  paper. 

Denote  by  <J>  (u,v)   the  generating  functions  of: 

A 

P  (x,y  1 1) ,  namely 

t   t 
Mu»v)  =  I   I   p  (x,y|t)  uV 
x=0  y=0 

Denote  $  as  the  generating  function  of  <J>  ,  then: 
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(1) 


oo 

*(u,v,e)  =  I    <j»t  et 


t=o 


£  0t  i?_       P(x,y|t)  uX  V8 


t    t 

I 

t=o   '  x=o  y=o 

«   [l-(qu  +  pv  +  ruv  +  s)  ] 

The  generating  functions  of  P(x,y|t)  can  be  expressed  as 


,-1 


(2) 


00        00/00 


*(U,V,6)   =  I  I       \l 

x=o  y=o  \t=max(x,y) 


P(x,y|t)9t  uX  vy 


*(u,v,e)  -  so  $(u,v,e) 


l  - 


l-s 


(qu  +  pr  +  ruv 


-1 


It  is  to  be  noticed  that  if  (x,y)il,  then  the  corresponding 
probabilities  p  and  q  in  the  above  expressions  do  not  give  the  required 
probabilities. 

III.   PROBABILITY  FUNCTIONS  FOR  THE  POINTS  ON  ABSORBING  BARRIERS 


Probabilities  P(x,y| t)  and  P(x,y| t)   can  be  computed  from 
formulas  (1)  and  (2).   In  order  _to  avoid  early  terminations  of  trials, 
these  probabilities  P(A)  and  P(B)  should  be  computed  and  evaluated  over 
sets  A  and  B  respectively,  as  illustrated  by  the  following  figure. 

(U) 


70% 


50% 


50%   70% 


These  probabilities  over  sets  A  and  B  should  be  very  close  to 
each  other  in  order  to  warrant  a  good  design  (i.e.,  trials  with  adequate 
data  points) . 
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If  P(A)  is  considerably  larger  than  P(B),  then  an  increase 
of  friendly  units  (or  decrease  of  the  threat  units)  in  the  original 
mix  is  necessary  to  guarantee  a  successful  trial.   If  P(A)  is  consid- 
erably smaller  than  P(B),  similar  policies  can  be  adopted. 


Denote  p(x,MJt)   the  absorption  probability  with  respect  to 
the  point   (x,M)eA  and  denote  p(N,y|t)   the  absorption  probability 
with  respect  to  the  point   (N,y)eB  at  the  end  of   t-th  (AT)   trial. 
The  expressions  of  these  absorption  probabilities  are  given  below: 


„,   1N,   x-k  M-k  k  t-x-M+k 
PfxMlO-  7       M(t-l)!  q E LJL 


(x-k)!  (M-k)!  k!  (t-x-M+k): 


(1) 


K(x,M,t) 

provided  0  <  x  <  N  and  M  >  1. 

«!/«.  m   N-k  y_lc  k  t-N-y+k  ... 

P(N  v|t)  =   Y        N(t-l).  q Vy      r  s y_ (2) 

K(N,Y|t)    (N"k):  (y_k):  kl  ^-N-y+k>: 

provided  0  <  y  <  M,    N  £  1. 

,      nN,   N-l-k  M-l-k  k  t+1+k-N-M    ,.. 

P(N  Mlt)  =  r   7  (t~1)'    q E E-S (3) 

r^,n|t;    r  k^   x       (N-l-k)!  (M-l-k)!  kl    (t-l+k-N+2-M)  I 

provided  N  >  1,   M  >  1. 

IV.  ABSORPTION  PROBABILITIES: 

Let  us  denote  P  (N,M)   the  probability  that  the  random  walk 
will  terminate  at  the  point   (a,M)  and  by  P  (N,M)   the  probability  that 

it  will  terminate  at  any  point  on  the  absorbing  barrier  A,  then  we  have: 

Po(N,M)  »  I  P(a,M|t) 

t=max(a,M) 

provided  0  <  a  <  N  and  M  >  1,  and 

N-l 
P  (N,M)  »  I  p  (N,M)   M  >  1. 

A       o=0    a 
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Denoting  by  t(ci,M)  the  number  of  updatings  to  be  performed  until  the 
particle  will  be  absorbed  at  point  x(a,M),  the  distribution  of   x(a,M) 
should  be  defined  by  the  conditional  probabilities  as  follows: 


d(t(<x  M)  =  t)  =  Pfo^l0 
pUU,n;   c;   p  (N)M) 

for  t  ■  max(a,M),  max(a,M)  +  1,  ... 

The  expected  number  of  updatings  is: 

t=max(a,M)    VaW»n> 


The  computation  (4)  is  pretty  combersome  for  the  above  quantity.   Yet 
E[t]   is  very  important  for  us  to  know.   If  it  turns  out  that  E[t], 
the  expected  number  of  updatings,  is  very  small,  we  will  have  insufficient 
data  during  the  trial.   If  E[t]  is  large,  then  the  selection  of   (AT) 
may  be  too  small  or  the  trial  may  last  too  long. 

The  following  paragraphs  will  be  devoted  to  a  discussion  how 
E[t]  can  be  used  in  computing  the  expected  duration  of  the  trial. 

If  (a,M)eA,  the  quantity: 

00 

I  t  P(o,M|t) 

t=max(a,M) 

a  M 

would  be  given  as  the  coefficient  of  the  term  u  v   in  the  expansion  of 

the  following  power  series: 

1  -'  (qu  +  s) 


1  -  (qu  +  pv  +  ruv  +  s) 


Denote  by  D.(N,M)   the  expected  duration  of  the  random  walk,  assuming  that 
it  will  terminate  at  any  point  on  the  absorbing  barrier  A.   Then  we  have: 

VN'M>  =  „  /VT  ^   I    E(T(a,M))  p  (N,M)   . 


A' 


VN'M)  a=0 


a 
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V.   APPLICATION 

Consider  the  rectangle  illustrated  in  Figure  3.   We  wish  to 
determine  the  mixes  of  forces,  i.e.,  to  determine  N,M  such  that  with 
high  probability  the  particle  will  be  absorbed  in  the  absorption  barrier 
in  the  set  C*  (see  Figure  3)  where: 

M  N        i 

C*  =  {(x,y):   either  x  =  N,  -  ^  y  <  M  or  y  =  M,  ^  ^  <  Nj  . 

In  other  words,  one  can  try  to  find  various  values  of  N,M  such  that  the 
absorption  probability  has  reached  to  some  acceptable  level  (say  80%, 
i.e.   P(C)  =  80%).   Once  such  an  optimal  mix  has  been  found,  experiments 
can'be  repeated  with  a  change  of  test  conditions,  i.e.,  terrains  or 
probability  of  kill  tables,  in  order  to  obtain  data  which  in  turn  will 
be  used  to  detect  any  probability  changes.   As  the  combat  between  the 
friendly  force  and  the  threat  force  going  on  in  the  battle  zone,  one 
keeps  a  continuous  record  on  the  quantity  defined  as  the  exchange  ratio 
at  time  t.   Definition  of  Exchange  Ratio  at  Time   t   (denote  as  ER(t))  is 


ER(t)  = 


// 

of  threat  tank  loss  up  to  and  including  time  t 

total  //  of  threat  tanks 

// 

of  friendly  anti-tank  systems  loss  up  to  and 
include  time  t 

total  //  of  friendly  anti-tank  systems 
(0,M) 


C       (N,M) 


p*     - 

i  y-  x  a  / 

— ^ 

* 

y 

jS 

(0,0) 


(N,0) 


Figure  3 


The  above  definition  is  used  for  "exchange  ratio  at  time  t"   if  data 
source  and  instruments  are  at  high  performance  level.   A  typical  graph 
of  this  exchange  ratio  ER(t)  is  given  in  Figure  4. 
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Figure  4. 


(0,0) 


(N,M) 


.  In  order  to  control  the  balance  of  force  structures  and 
mixes,  It  is  ideal  to  have  some  bounds,  say  a  and  b  such  that 
the  probability  that  ER(t)   is  bounded  in  between  a  and  b   for 
all  t  >  0  is  very  high,  i.e., 


P  {a  <  ER(t)  <  b;  0  <  t}  =  .95 

It  is  clear  that  a  =  tan  8,   and  b  =  tan  6-  where  9,   and  89 

are  shown  in  Figure  4. 

Using  a  two-dimensional  random  walk  model  [2. 4],  we  have  illus- 
trated above  how  to  determine  the  best  force  structures  such  that  the 
exchange  ratio  is  close  to  one  for  most  of  the  times  during  a  trial.   In 
other  words,  the  probability  that  the  sample  path  falls  into  the  shaded 
region  (see  Figure  5)  is  maximized  by  an  optimal  choice  of  force  structures 


(0,0) 


Figure  5 


(1,D 


Due  to  the  fact  that  instruments  are  most  likely  working  at 
low  performance  levels,  we  provide  an  alternative  definition  for  "exchange 


-  10  - 


ratio"  using  indirect  measurements.   This  proposed  definition  of  "exchange 
ratio"  is  the  ratio  of  expected  number  of  kills  normalized  by  the  initial 
force  units. 

The  above  definition  of  the  exchange  ratio  is  closer  to  the 
measure  of  effectiveness  rather  than  the  exchange  of  kill  ratio  but  due 
to  the  instrument  limitations  this  may  be  a  good  alternative  measure  of 
exchange  ratio. 

Another  way  of  defining  the  same  quantity  is  by  looking  at  the 
remaining  kill  potential  or  fighting  value  of  the  forces.   This  remaining 
fighting  force  would  be  compared  with  their  respective  fighting  potentials 
relative  to  their  respective  initial  fighting  potential.   The  degradation 
of  their  respective  potentials  over  time  will  be  plotted  on  a  graph  sheet 
for  each  trial. 

The  literature  indicated  that  discrete  state  stochastic  model 
formulations  of  combat  have  been  difficult  to  solve  even  when  the  process 
is  considered  to  be  Poisson  (Lanchester  type)  with  stationary  transistion 
mechanisms  JlJ  . 

The  technical  portion  of  our  discussion  seems  new  and  provides 
a  discrete  state  two-dimensional  stochastic  model  of  combat  which  can  be 
extended  into  a  complete  stochastic  model  for  field  experimentation  such 
as  the  Tactical  Effectiveness  Testing  of  Anti  tank  Missiles  [3J. 
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